Pressure and linear heat capacity in the superconducting state of thoriated UBei 3 



R.J. Zieve 1 , R. Duke 1 and J.L. Smith 2 

1 Physics Department, University of California at Davis 
2 Los Alamos National Laboratory, Division of Materials Science and Technology 

Even well below T c , the heavy-fermion superconductor (U,Th)Bei3 has a large linear term in its specific 
heat. We show that under uniaxial pressure, the linear heat capacity increases in magnitude by more than a 
factor of two. The change is reversible and suggests that the linear term is an intrinsic property of the material. 
In addition, we find no evidence of hysteresis or of latent heat in the low-temperature and low-pressure portion 
of the phase diagram, showing that all transitions in this region are second order. 
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For the past two decades, heavy-fermion superconductors 
have revealed a variety of unusual behaviors that hint at un- 
conventional superconductivity |1]. Many low-temperature 
properties such as specific heat and NMR relaxation rates have 
power law rather than exponential temperature dependences. 
Two compounds, UPt3 and Ui-ajTh^Beis for x between 0.02 
and 0.04, each have two transitions leading to distinct su- 
perconducting phases. Further phase transitions appear with 
pressure or applied magnetic field. Yet no experiment has 
emerged that conclusively identifies the order parameters of 
the phases or even their symmetry. 

One potential clue to heavy-fermion order parameters is 
the significant linear term in the thermal properties within 
the superconducting state. In theory linear heat capacity 
and thermal conductivity are normal-state phenomena, which 
should disappear once the superconducting energy gap alters 
the excitation spectrum. Yet specific heat C(T) — ^ S T + 
C non -ii ne ar {T) in the superconducting phase is a persistent 
feature in heavy-fermion superconductors. The coefficient 7 S 
can reach over 50% of j n , the normal-state value of C(T)/T 
just above the transition. First seen in UPt3 [2], a large 7 S 
and an analogous linear term in thermal conductivity are also 
found in CeCoIn 5 0, URu 2 Si 2 fl H, and UPd 2 Al 3 ji H, 
among others. 

The linear term is sometimes viewed as stemming entirely 
or in part from imperfect samples. One simple explanation 
would be that some fraction of the material remains normal. 
Another heavily discussed alternative, resonant impurity scat- 
tering, combines impurity effects with intrinsic features of 
the order parameter |§l 1 1 Oil - For a d-wave superconduc- 
tor, a very small impurity concentration could create finite un- 
gapped regions on the Fermi surface and a constant density 
of states. The strength, phase shift, and anisotropy of the im- 
purity scattering and the impurity density all factor into the 
normal-like behavior. Yet other scenarios treat the linear term 
as an intrinsic property. These include involvement of only 
part of the Fermi surface in superconductivity, and the more 
exotic odd-frequency pairing fflTlll . 

In support of the importance of impurities, the magnitude of 
7s varies significantly from sample to sample, ranging from 
0.127„ to 0.627„ in UPt3. Furthermore, 7 S generally de- 
creases as the superconducting transition temperature of the 
sample increases. Since a higher T c often indicates a better- 
quality sample, the correlation does suggest that j s arises at 



hole for 
sample access 



sample 




fill line 



helium bellows 



piezo sensor 
NbTi spacers 



FIG. 1 : Helium bellows setup for measuring specific heat under uni- 
axial pressure. 



least in part from sample problems |2j]. Also, j s in pure 
UBei3 is small or zero, but doping with either thorium fl2ll 
or boron II 3lll411 increases both 7 S and the likelihood of sam- 
ple inhomogeneities. 

Any theory of the linear terms in the superconducting spe- 
cific heat and thermal conductivity must also address the ab- 
sence of a linear contribution to the NMR spin-lattice relax- 
ation rate. If the large j s comes from either a normal portion 
of the sample or a finite region of the Fermi surface with no 
gap, that same source should lead to a linear Korringa relax- 
ation, as in the normal state. However, NMR measurements 
on a variety of heavy-fermion superconductors find only a cu- 
bic temperature dependence down to temperatures well below 
where linear terms in heat capacity and thermal conductivity 
become significant 1 15lll6tll7ll . The cubic dependence would 
be expected from line nodes in the gap. Reconciling the NMR 
and thermodynamic data would certainly be a step towards 
understanding heavy-fermion superconductivity. 

Our present work shows a large and reversible change in 
7s with pressure in Uo.98Tho.o2Bei 3 . Explaining the pressure 
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FIG. 2: Specific heat of Uo.9sTho.o2Bei3 as a function of tempera- 
ture, for several applied uniaxial pressures. Inset enlarges part of the 
low-temperature region, showing that the curves for different pres- 
sures are nearly parallel. 



FIG. 3: Linear coefficient of specific heat in superconducting phase, 
7s, as a function of pressure. Inset: Best fit parameters for low 
temperature tail: the exponent n from a 3-parameter fit C(T) = 
-jsT + AT" (top) and the pref actor A from a 2-parameter fit C(T) = 
7 S T + AT Z (bottom). 



dependence, which occurs without a change in the impurity 
concentration, will further restrict theoretical treatments. 

We use a pressure cell activated by a helium bellows. The 
setup is mounted at the mixing chamber of a KelvinOx 100 di- 
lution refrigerator, and we can change pressure while keeping 
the sample temperature below 300 mK. Our cell, illustrated 
schematically in Figure^ is modelled after the cell described 
by Pfleiderer et al. Hal . The expanding bellows presses on 
a column including the sample and a piezoelectric crystal to 
measure pressure changes. The small cross-sectional area of 
the sample amplifies the pressure within the bellows; by the 
time the helium solidifies at 25 bar we reach a uniaxial pres- 
sure of 7.8 kbar at the sample. The uniaxial technique is nec- 
essary for changing pressure at low temperature. To avoid 
symmetry -breaking effects from the uniaxial pressure, we use 
a polycrystalline sample in the experiment. 

We use a transient pulse method for heat capacity measure- 
ments. Our heater is a 50:50 AuCr thin film, our thermome- 
ter a RUO2 film. Pieces of NbTi on each side of the sample 
provide a thermal link, with a time constant of order 8 sec- 
onds between the sample and the rest of the bellows. As a 
conventional superconductor well below T c , the NbTi itself 
contributes negligibly to the measured heat capacity. 

One difficulty with the measurements is a significant time 
constant between the helium bellows and the dilution refriger- 
ator, of order 3 minutes near 500 mK and increasing to 10 
minutes at 200 mK. Waiting for the bellows to equilibrate 
completely with the cryostat at each temperature takes pro- 
hibitively long. Instead, we measure the relaxation time of the 
bellows throughout our temperature range and verify that it is 
independent of pressure. We then account for a slowly chang- 



ing bellows temperature in our fits of the temperature decays 
after each heat pulse. 

In Figure [2] we show C/T as a function of temperature 
for five different pressures. At the lowest pressure there are 
two transitions, centered at 450 mK and 550 mK, with a small 
shoulder between them. As pressure increases, the two transi- 
tions decrease in temperature and merge into one, the shoulder 
disappearing. The amplitude of the peak decreases, while that 
of the low-temperature tail increases. The normal-state heat 
capacity does not change with pressure. All this agrees with 
the previous uniaxial pressure experiment fl9ll . 

The lowest-pressure data follows the power-law form pre- 
viously observed for (U,Th)Bei 3 , C(T) = j s T + AT 3 with 
7 S = 0.64 J/mol K at P = 0.03 kbar. Previous low- 
temperature heat capacity measurements for different Th con- 
centrations found C(T) = 7s T + AT n with a best-fit 
exponent n near 4, an unphysical value, for Th concentra- 
tions of 2.2% and above. Our sample, with a slightly lower Th 
concentration, retains the T 3 behavior. Indeed, this function 
proves to fit our data well for all pressures. Three-parameter 
fits of C(T) to the form j s T + AT n give an exponent n w 3 
at all pressures, as shown in Figure [3] With this in mind, we 
fix n = 3 and carry out two-parameter fits. We find a steady 
increase in 7 S with pressure, while A decreases more slowly; 
these quantities are shown in Figure[3] 

As a further check, we extrapolate the specific heat curves 
to T = according to the above fits. We then integrate to find 
the entropy S by S(T) = f Q T dT^P-. FigureHcompares the 
entropy for the lowest and highest pressure curves of Figure|2] 
By 700 mK, safely in the normal state, the total entropy varies 
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FIG. 4: Entropy from to 700 mK for 0.03 kbar (solid) and 5.58 
kbar (dashed), the lowest and highest pressures in Figure [2] Inset: 
Specific heat jump at the transition, for single-transition regime. 



by only a few percent among pressures. The extra entropy 
under the C/T peak at low pressures offsets the extra entropy 
under the low-temperature tail at high pressures. 

With this confirmation, we return to the j s . The magnitude 
of the change with pressure is striking, more than a factor of 
two. Furthermore, the change is completely reversible, even 
without raising the temperature above 300 mK. In fact, ex- 
ploring the reversibility of the low-temperature heat capacity 
originally motivated our measurements. The phase diagram of 
Ui_a;Th x Bei3 includes boundaries within the superconduct- 
ing regime as functions of concentration and pressure, as well 
as temperature. For Th concentrations with two transitions, 
/iSR measurements finds weak local magnetic order below 
the lower-temperature transition 1 20|, 12111 . No local order ap- 
pears for Th concentrations with a single transition, suggest- 
ing a phase boundary near x = 0.02 between the single phase 
at x < 0.02 and the lower-temperature phase for x > 0.02. 
Since increasing pressure acts much like decreasing Th con- 
centration I 22II23I1 . pressure measurements can cross an anal- 
ogous phase boundary and explore its thermodynamic proper- 
ties. 

The earlier heat capacity measurements under pressure lll9ll 
drove the sample around rather than across the phase bound- 
ary: pressure was always changed at room temperature, with 
the sample then cooled into one phase or the other. Our bel- 
lows pressure cell allows us to change pressure while cold, 
thereby crossing the transition directly. We also use a sam- 
ple with Th concentration closer to the transition and extend 
the measurements to lower temperatures by using a dilution 
refrigerator rather than the pumped 3 He cryostat of the earlier 
work. 

We use several paths through the pressure-temperature 



space to search for hysteresis. In one case, we change pres- 
sure, keeping temperature below 300 mK. We then measure 
C(T) from low temperature to above T c . After the sample 
has warmed above T c , we cool and repeat the specific heat 
measurements from our lowest temperatures. The two spe- 
cific heat curves agree to better than 0.5% for both increasing 
and decreasing pressure, up to 5.5 kbar. Changing pressure 
at fixed temperature is less reliable, both because it is more 
difficult to return exactly to the original pressure and because 
the temperature always changes slightly during the pressure 
change, but again there is no evidence of hysteresis in the spe- 
cific heat. 

We also find no evidence of latent heat when we monitor the 
sample temperature while increasing or decreasing pressure. 
The temperature rises during any pressure change. The effect 
is more noticeable upon increasing pressure, but the difference 
is completely explained by the heat load to the cryostat from 
adding additional room-temperature helium and by the work 
done in compressing the sample. We never find a tempera- 
ture reduction on changing pressure. We conclude that any 
transition with pressure in this region is second-order. 

Although pressure might introduce additional defects into 
a sample, perhaps even defects that could substantially alter 
the heat capacity, such an explanation for the change in 7^ 
also demands, implausibly, that the defects anneal away at low 
temperature. The large reversible effect on 7 S suggests that in 
fact the 7 S is tied intimately to the mechanism of supercon- 
ductivity itself. 

Note that at our highest pressures 7 S exceeds the normal- 
state C/T just above T c . This confirms other evidence of 
non-Fermi liquid behavior in UBei3. An entropy deficit in 
normal UBei3 relative to the superconducting phase has long 
been known, suggesting that j n (T) increases substantially be- 
low T c . Suppressing T c with a magnetic field bears this out, 
with C/T increasing steadily toward lower temperatures. The 
same effect, with an even stronger increase in j n (T), appears 
for Th-doped UBei3. Heat capacity in a 3 T field shows j n (T) 
of 1400 mJ/mole K 2 at 0.42 K. To match the measured super- 
conducting entropy it must rise to 2300 mJ/mole K 2 at T = 0, 
a faster than linear increase in 7„ (T) itself El. 

The behavior of j a emphasizes that increasing pressure and 
decreasing Th concentration have analogous but not identi- 
cal effects. As shown previously, the topology of the phase 
diagram appears to be similar for the two variables, but the 
temperature dependence of the transitions is not. On decreas- 
ing Th concentration below 2%, T c rises. While pressure also 
merges the transitions, T c decreases monotonically. Similarly, 
7 S generally decreases with decreasing Th concentration; it is 
more an order of magnitude smaller in pure UBei3 than for 
2% Th doping. Yet 7^ increases with increasing pressure. 

This suggests that changes in 7^ come from quantitative 
rather than qualitative changes in the order parameter. For 
a BCS superconductor, the specific heat jump at the transi- 
tion, AC, and the magnitude of the energy gap, A(0), satisfy 
AC = 1.437„T C and A(0) = 1.76kT c . Although these sim- 
ple proportionalities fail for a non-s-wave order parameter or 
strong coupling, the discontinuity is still related to A(0). We 
fit our specific heat data with a single sharp transition, with en- 
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tropy conserved between our data and the fit. The inset of Fig- 
ure@]shows the size of the jump, at pressures high enough that 
a single transition provides a good fit to the data. The jump 
decreases with increasing pressure, although not as rapidly as 
7s increases. 

Other heavy-fermion materials also show substantial 
changes in 7 S with pressure. In UPd2Al3, with supercon- 
ducting T c near 1.5K and antiferromagnetic Tjy ~ 18K, j s 
increases over 50% at 10.8 kbar of hydrostatic pressure |6]. 
In this case the heat capacity just above T c , deep in the an- 
tiferromagnetic phase, increases by a comparable amount. A 
possible explanation is that separate electron subsystems are 
responsible for the magnetic and superconducting behaviors, 
with the pressure dependence arising only from the electrons 
responsible for the magnetism. 

Another example is CeRhlns (25il at the boundary between 
antiferromagnet and superconductor. In this case, the super- 
conducting 7 S vanishes at high pressures, but rises steadily as 
pressure decreases from 21 to 15 kbar. At 15 kbar, at the an- 
tiferromagnetic transition, 7 S has reached its value within the 
antiferromagnet. The authors interpret the large linear term 
as a consequence of finite regions of ungapped Fermi surface 



ll25ll that come from an increase in anisotropic impurity scat- 
tering near the antiferromagnetic transition [ 1011 . 

Thus a large and changing 7 S appears not only in our su- 
perconducting system but also in a superconducting antifer- 
romagnet and at a superconductor/antiferromagnet transition. 
Whether there is any further connection among these systems, 
such as proximity to an unrealized antiferromagnetic transi- 
tion in (U,Th)Bei3, remains to be seen. 

In summary, we observe a large increase with pressure in 
7 S , the linear coefficient of the superconducting specific heat. 
We also find that all pressure-dependent behavior is reversible, 
indicating that all phase transitions in the region are second or- 
der. The change in 7 S with no change in the impurity density 
appears inconsistent with some proposed explanations for the 
origin of the linear term, including resonant impurity scatter- 
ing. Whether or not a non-zero 7 S is itself an intrinsic property 
of (U,Th)Bei3, its strong variation within a single sample is 
likely intrinsic and may prove a useful signature of the nature 
of superconductivity in the material. 
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